Charged surface in salty water with multivalent ions: Giant inversion of charge. 
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Screening of a strongly charged macroion by oppositely charged colloidal particles, micelles, or 
short polyelectrolytes is considered. Due to strong lateral repulsion such multivalent counterions 
form a strongly correlated liquid at the surface of the macroion. This liquid provides correlation 
induced attraction of multivalent counterions to the macroion surface. As a result even a moderate 
concentration of multivalent counterions in the solution inverts the sign of the net macroion charge. 
We show that at high concentration of monovalent salt the absolute value of inverted charge can be 
larger than the bare one. This giant inversion of charge can be observed in electrophoresis. 
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Charge inversion is a phenomenon in which a charged 
particle (a macroion) strongly binds so many counterions 
in a water solution, that its net charge changes sign. As 
shown below the binding energy of counterion with large 
charge Z is larger than ksT, so that this net charge is 
easily observable; for instance, it is the net charge that 
determines linear transport properties, such as particle 
drift in a weak fiflld electrophoresis. Charge inversion 
has been observed^ in polyelectrolyte-micelle system and 
is possible for a variety of other systems, ranging from 
solid surface of mica or lipid membranes, to DNA or 
actin. 

Charge inversion is of special interest for delivery of 
genes to the living cell for the purpose of gene therapy. 
The problem is that both bare DNA and a cell surface 
are negatively charged and repel each other. The goal is 
to screen DNA in such a way that the resulting complex 
is positivecl. 

Theoretically, charge inversion can be also thought 
of as an over-screening. Indeed, the simplest screening 
atmosphere, familiar from linear Dcbye-Hiickel theory, 
compensates at any finite distance only a part of the 
macroion charge. It can be proven that this property 
holds also in non-linear Poisson-Boltzmann (PB) theory. 
The statement that the net charge preserves sign of the 
bare charge agrees with the common sense. One can 
think that this statement is even more, universal than re- 
sults of PB equation. It was showntTu, however, that 
this presumption of common sense fails for screening 
by Z-valent counterions (Z-ions), such as charged col- 
loidal particles, micelles, or short polyelectrolytes, be- 
cause there are strong lateral correlations between them 
when they are bound to the surface of a macroion. These 
correlations are beyond the mean field PB theory, and 
charge inversion is their most spectacular manifestation. 

Charge inversion has attracted a significant attention 
in the last couple of yearsH. Our goal in the present pa- 
per is to provide a simple physical explanation of charge 
inversion and to show that in the most practical case, 
when both Z-ions and monovalent salt, such as NaCl, are 
present, not only charge sign may flip, but the inverted 
charge can become even larger in absolute value than the 
bare charge, thus giving rise to giant charge inversion. 



Let us demonstrate the role of lateral correlations be- 
tween Z-ions for a primitive toy model. Imagine a hard- 
core sphere with radius b and with negative charge Q 
screened by two spherical positive Z-ions with radius a. 
One can see that if Coulomb repulsion between Z-ions 
is much larger than k B T they are situated on opposite 
sides of the negative sphere (Fig. la). If Ze < 2\Q\ each 
Z-ion is bound, because the energy required to remove 
it to infinity \Q\Ze/(a + b) — Z 2 e 2 /2(a + b) is positive. 
Thus, the charge of the whole complex Q + 2Ze can be 
positive and as large as 3|Q|. This example demonstrates 
the possibility of an almost 300% charge inversion. It is 
obvious that this charge inversion is a result of the corre- 
lation between Z-ions which avoid each other and reside 
on opposite sides of the negative charge. On the other 
hand, description of screening of the central sphere in PB 
approximation smears the positive charge, as shown on 
Fig. lb and does not lead to the charge inversion. Indeed, 
in this case charge accumulates in spherically symmetric 
screening atmosphere only until the point of neutrality 
at which electric field reverses its sign and attraction is 
replaced by repulsion. 
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FIG. 1. a) A toy model of charge inversion, b) PB approx- 
imation does not lead to charge inversion. 

In this paper we consider screening of a macroion sur- 
face with negative immobile surface charge density — a 
by finite concentration of positive Z-ions, neutralizing 
amount of monovalent coions, and a large concentration 
Ni of a monovalent salt. This is more practical problem 
than one considered in Ref. where monovalent salt 
was absent. Correspondingly, we assume that all inter- 
actions are screened with Debye-Hiickel screening length 
r s = (87r2 B 7Vi) _1/2 , w here l B = e 2 /{Dk B T) is the Bjer- 
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rum length, e is the charge of a proton, D ~ 80 is the 
dielectric constant of water. 

We begin with the simplest macroion which is a thin 
charged sheet immersed in water solution (Fig. 2a). 
Later we examine more realistic macroion which is a thick 
insulator charged at the surface (Fig. 2b). 
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FIG. 2. Models studied in this paper. Z-ions are shown by 
full circles, a) Charged plane immersed in water. b)Surface of 
a large macroion. Image charges are shown by broken circles. 

Assume that the plane with the charge density — a is 
covered by Z-ions with two-dimensional concentration n. 
Integrating out all monovalent ions, or, equivalently, con- 
sidering all interactions screened at the distance r s , we 
can write down the free energy per unit area in the form 



F = 7T(T 2 r s /L> - 2nar s Zen/D + F zz + F, 



id : 



(1) 



where the four terms are responsible, respectively, for the 
self interaction of the charged plane, for the interaction 
between Z-ions and the plane, for the interaction between 
Z-ions and for the entropy of ideal two-dimensional gas 
Z-ions. 

Our goal is to calculate the net charge density of the 
plane 



a* = —a + Zen. 
Using Eq. (||) one can rewrite Eq. (Q) as 
F = ir{a*) 2 r s /D + F OC p, 



(2) 
(3) 



where Focp — F c + Fid is the free energy of the same 
system of Z-ions residing on a neutralizing background 
with surface charge density —Zen, which is convention- 
ally referred to as one component plazma (OCP), and 
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-n{Zen) 2 r s /D + F, 
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(4) 



is the correlation part of Focp ■ This transformation can 
be simply interpreted as the addition of uniform charge 
densities —a* and a* to the plane. The first addition 
makes a neutral OCP on the plane. The second plane of 
charge creates two plane capacitors with negative charges 
on both sides of the plane which screen inverted charge 



of the plane at the distance r s . The first term of Eq. (g) 
is nothing but the energy of these two capacitors. There 
is no cross term in energy between the OCP and the ca- 
pacitors because each plane capacitor creates a constant 
potential, ip(0) = 2na*r s /D, at the neutral OCP. 

Using Eq. (Q), the electrochemical potential of Z-ions 
at the plane can be written as fj, = Zeip(0) + fiid + /i C! 
where \iid and u c = dF c /dn are the ideal and the cor- 
relation parts of the chemical potential of OCP. In equi- 
librium, /i is equal to the chemical potential, /i& of the 
bulk solution, because in the bulk electrostatic potential 
tp = 0. Using Eq. (||), we have: 



2na*r s Ze/D = -fj, c + (p b - fi td ) 



(5) 



As we show below, in most practical cases the correla- 
tion effect is rather strong, so that [i c is negative and 
|Mc| 3> ksT. This means that for large enough con- 
centration of Z-ions in the bulk and at the surface, n, 
both bulk chemical potential fib and ideal part of surface 
chemical potential /i^ should be neglected compared to 
fi c . Furthermore, strong correlations imply that at least 
short range order of Z-ions on the surface should be sim- 
ilar to that of triangular Wigner crystal (WC) since it 
delivers the lowest energy to OCP. Therefore, 
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We see now that the net charge density a* is positive. 
This proves inversion of the bare charge density — a . Eq. 
(||) has a very simple meaning: \nwc\/Ze is the "corre- 
lation" voltage which charges two above mentioned par- 
allel capacitors with thickness r s and total capacitance 
per unit area D/(2-Kr s ). 

To calculate the "correlation" voltage \nwc\/Ze, we 
start from the case of weak screening when r s is larger 
than the average distance between Z-ions. In this case, 
screening does not affect thermodynamic properties of 
WC. The energy per Z-ion e(n) of such Coulomb WC at 
T = can be estimated as an interaction energy of a Z- 
ion with its Wigner-Seitz cell, because interaction energy 
of neigboring neutral Wigner-Seitz cells is very small. 
This gives e(n) = —Z 2 e 2 /RD, where R = (to)" 1 ^ is the 
radius of a Wigner-Seitz celLdwe approximate hexagon 
by a disc). More accurately!] e(n) = — \.\Z 2 e 2 j RD = 
-lMn 1 / 2 Z 2 e 2 /D. One can discuss the role of a finite 
temperature on WC in terms of the inverse dimension- 
less temperature T = Z 2 e 2 / (RDksT). We are inter- 
ested in the case of large T. For example, at a typical 
Zen = a = 1.0 e/nm 2 and at room temperature, T = 10 
even for Z = 4. Wigner crystal meltafl at T = 130, so 
that for r < 130 we deal with a strongly correlated liq- 
uid. Numerical calculations, however, confirm that at 
r ^> 1 thermodynamic properties of strongly correlated 
liquid are close to that of WCh. Therefore, for estimates 
of fi c we can still write that F c — ne(n) and use 



d (ne(n)) 
dn 



= -imVk B T = -1.65 



Z 2 e 2 
RD 



(7) 
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We see now that indeed fJ,wc is negative and |jtiwc| ^> 
k B T, so that Eq. (^|) is justified. Substituting Eq. (0) 
into Eq. (§), we get a* = 0.83Ze/(irr s R). At r s > R, 
charge density a* <C <r, and Zen ~ a, one can replace R 
by i? = (cTTr/Ze)- 1 / 2 . This gives 

a* /a = 0.83(R /r s ) = 0.83( 1/2 , (C « 1) (8) 

where £ = Ze/irar 2 is a dimensionless charge of a Z-ion. 
Thus, at r s ^ i? or £ <C 1, inverted charge density grows 
with decreasing r s . Extrapolating to r s = 2Rq where 
screening starts to substantially modify the interaction 
between Z-ions we obtain a* = OAa. 

Now we switch to the case of strong screening, r s <C R, 
or £ ^> 1. It seems that in this case a* should decrease 
with decreasing r s , because screening reduces the energy 
of WC and leads to its melting. In fact, this is what 
eventually happens. However, there is a range of r s <C R 
where the energy of WC is still large. In this range, as r s 
decreases, the repulsion between Z-ions becomes weaker, 
what in turn makes it easier to pack more of them on the 
plane. Therefore, a* continues to grow with decreasing 
r s - 

At r s <C R one is still able to estimate thermodynamic 
properties of OCP from the model of a triangular WC. 
Keeping only interactions with the 6 nearest neighbors in 
Eq. (jij), we can write the correlation part of free energy 
of screened WC per unit area as 

■nr s (Zen) 2 (Ze) 2 , . , . 

Fc = ^^ + 3n^-exp(-A/r s ), (9) 

where A = (2/y / 3) 1 / 2 n^ 1 / 2 is the lattice constant 
of this WC. Calculating the chemical potential of Z- 
ions at the plane, [iwc — dF c /dn and substituting 
it into Eq. (|) one finds that A ~ r s ln(3C/4), R ~ 
(27r/v / 3) 1/2 r s ln(3C/4) and 

— = F 2n 2 C - 1, (C»i). (io) 

o V3 ln 2 (3C/4) 

Alternatively, one can derive Eq. ([l^) by direct mini- 
mization of Eq. (|l]) with respect of n. In this way, one 
does not need a capacitor interpretation which is not as 
transparent in this case as for r s ^S> R- 

Thus, at r s <C R, or ( 3> 1 the distance R decreases 
and inverted charge continues to grow with decreasing 
r s . This result could be anticipated for the toy model 
of Fig. la if Coulomb interaction betwen the spheres is 
replaced by a strongly screened one. Screening obviously 
affects repulsion between positive spheres stronger than 
their attraction to the negative one and, therefore, makes 
maximum allowed charges Ze larger. 

Above we studied analytically two extremes, r s R 
and r s <ti R. In the case of arbitrary r s we can find 
a* numerically. For this purpose we calculate /iwc from 
Eq. (|) and substitute it in Eq. (g). This gives 



where the sum is taken over all vectors of WC lattice 
and can be evaluated numerically. Then one can find the 
equilibrium n for any given values of The resulting 
ratio a* I a is plotted by the solid curve in Fig. 3. 
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FIG. 3. The ratio a* /a as a function of the charge £. The 
solid curve is calculated for a charged plane by a numerical 
solution to eq. the dashed curve is the large r s limit, 

eq. (ph. The • points are calculated for the screening of the 
surface of the semispace with dielectric constant much smaller 
than 80. In this case image charges (Fig. 2b) are taken into 
account. 



Since the value of a* represents the main result of our 
work, its subtle physical meaning should be clearly un- 
derstood. Indeed, the entire system, macroion plus over- 
charging Z-ions, is of course neutralized by the monova- 
lent salt. One can ask then, what is the meaning of charge 
inversion? The answer is simple for r s 3> R, when charge 
a* is well separated in space from the oppositely charged 
atmosphere of monovalent salt (which leads to the in- 
terpretation based on two capacitors, see above). When 
r s <C R there is no such obvious spatial separation. Nev- 
ertheless, a* can be observed, because Z-ions are bound 
with energies well above fc^T while small ions are only 
weakly bound. First, the number of bound Z-ions can 
be counted using, e.g., the atomic force microscopy. Pos- 
itive a* means " over-population" : there are more bound 
Z-ions than neutrality condition implies. Second, it is 
a* that determines the mobility of macroion in the weak 
field electrophoresis experiments. 

The results discussed so far were derived for the 
charged plane which is immersed in water and screened 
on both sides by Z-ions and monovalent salt (Fig. 2a). 
In reality charged plane is typically a surface of a rather 
thick membrane whose (organic, fatty) material is a di- 
electric with permeability much less than that of water. 
In this case, image charges which have the same sign as 
Z-ions must be taken into account (Fig. 2b). We have 
analyzed this situation in details, which will be reported 
elsewhere. The main result turns out to be very simple: 
while image charges repel Z-ions and drive the entire 
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Wigner crystal somewhat away from the surface, their 
major effect is that in this case only one capacitor must 
be charged (on the water side of the surface). Accord- 
ingly, the ratio a* /a is reduced by a factor very close to 
2 compared to the case of two-sided plane (Fig. 3). 

We are prepared to address now the question of max- 
imal possible charge inversion. How far can a macroion 
be overcharged, and what should one do to achieve that? 
Figure 3 and equation (9) suggest that the ratio a* jo 
continues to grow with growing £. However, the possibil- 
ities to increase C are limited along with the assumptions 
of the presented theory. Indeed, there are two ways to 
increase C — Ze/airr 2 , namely to choose surface with 
small a and ions with large Z. The former way is re- 
stricted because Z-ions remain strongly bound to the 
surface only as long as I/^vkcI — 2irr s aZe/D 3> fc^T 
or £ < 2Z 2 ls/r s . Therefore, the latter way, which is 
to increase Z, is really the most important. It is, how- 
ever, also restricted, because atJarge Z , monovalent ions 
start to condense on the Z-ior£9. Assuming Z-ions are 
spheres of the radius a, their effective net charge at large 
Z can be written as Z e g = (a/Is) 2 In (ZZsr s /a 2 ), yield- 
ing C < 8 (a 2 /Zsr s ) [in (Z^r s /a 2 )] 2 . Since this estimate 
was derived under the assumption that r s > a, the largest 
a we can choose is a = r s . For r s = a = 10 A charge £ 
may be as high as about 10, so that the ratio a* jo can 
exceed 100%. 

Since charge inversion grows with increasing a we are 
tempted to explore the case a > r s . To address this 
situation, our theory needs a couple of modifications. 
Specifically, in the first term of Eq. (^) we must take 
into account the fact that only a part of Z-ion interacts 
with the surface, namely the segment which is within the 
distance r s from the surface. One should also take into 
account that strong screening increases Z e ff. Assuming 
Z-ion is a sphere, this modifies upper bound for £ by 
a factor a/r s and thus it makes charge inversion even 
larger. We do not discuss this regime in details, because 
it is highly non-universal, dependent on the shape and 
charge distribution of the Z-ions, plane roughness, etc. 

Meanwhile, there is much more powerful way to in- 
crease charge inversion. Suppose we take Z-ions with the 
shape of long rigid rods. Such a situation is very prac- 
tical, since it corresponds to the screening of charged 
surface by rigid polyelectrolytes, such as DNA double 
helixO. In this case, correlation between Z-ions leads to 
parallel, nematic-like ordering of rods on the surface. In 
other words, WC in this case is one-dimensional, perpen- 
dicular to rods. Chemical potential |/iiyc| in this case is 
about the interaction energy of one rod with the stripe 
of the surface charge, which plays the role of the Wigner- 
Seitz cell. Importantly, this energy, along with the effec- 
tive net charge, Z e s , are proportional to the rod length L 
and thus can be very large. Rods can be strongly bound, 
with chemical potential much exceeding ksT, even at 
very small a. This holds even in spite of the Onsager- 
Manning condensation!! 2 ] of monovalent ions on the rods: 



for instance, at A > r s > a one has Z e ff = Lrj c /e , where 
A and a are, respectively, the distance between rods in 
WC and radius of the rod (double helix), r/ c — ksT/e. 
As a result the ratio a* /a grows with decreasing r s as 
a* I a ~ (i] c /2r s a) In (rj c /2TTr s a). At r s ~ a and small 
enough a this ratio can be much larger than one. This 
phenomenon can be called giant charge inversion. 

Giant charge inversion can be also achieved if DNA 
screens a positively charged wide cylinder with the radius 
greater or about the DNA double helix persistence length 
(500 A). In this case DNA spirals around the cylinder, 
once again with WC type strong correlations between 
subsequent turns. We leave open the possibility to spec- 
ulate on the relevance of this model system to the, fact 
that DNA overcharges a nucleosome by about 20%fl. 

To conclude, we have presented simple physical argu- 
ments explaining the nature and limitations of charge 
inversion in the system, where no interactions are oper- 
ational except for Coulomb and short range hard core 
repulsion. Correlations between bound ions, which are 
strong for multivalent counterions with Z » 1, are the 
powerful source of charge inversion for purely electro- 
static system. We have shown that even spherical Z-ions 
adsorbed on a large plane macroion can lead to charge in- 
version larger than 100%, while for rod-like Z-ions charge 
inversion can reach gigantic proportions. 
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1 Y. Wang, K. Kimura, Q. Huang, P. L. Dubin, W. Jaeger, 
Macromolecules, 32 (21), 7128 (1999). 

2 P. L. Feigner, Sci. American, 276, 86 (1997). 

3 J. Ennis, S. Marcelja and R. Kjellander, Electrochim. Acta, 
41, 2115 (1996) 

4 V. I. Perel and B. I. Shklovskii, Physica A 274, 446 (1999). 

5 B. I. Shklovskii, Phys. Rev. E60, 5802 (1999). 

6 E. M. Mateescu, C. Jeppersen and P. Pincus, Europhys. 
Lett. 46, 454 (1999); S. Y. Park, R. F. Bminsma, and W. 
M. Gelbart, Europhys. Lett. 46, 493 (1999); J. F. Joanny, 
Europ. J. Phys. B 9 117 (1999). 

7 L. Bonsall, A. A. Maradudin, Phys. Rev. B15, 1959 (1977). 

8 R. C. Gann, S. Chakravarty, and G. V. Chester, Phys. Rev. 
B 20, 326 (1979). 

9 H. Totsuji, Phys. Rev. A 17, 399 (1978) 

10 M. Guerom, G. Weisbuch, Biopolimers, 19, 353 (1980); S. 
Alexander, P. M. Chaikin, P. Grant, G. J. Morales, P. Pin- 
cus, and D. Hone, J. Chem. Phys. 80, 5776 (1984); S. A. 
Safran, P. A. Pincus, M. E. Cates, F. C. MacKintosh, J. 
Phys. (France) 51, 503 (1990). 

11 Ye Fang, Jie Yang, J. Phys. Chem. B 101, 441 (1997). 

12 G. S. Manning, J. Chem. Phys. 51, 924 (1969). 



4 



